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Abstract 

It is shown that after reflection from a lossy medium the s- or p- 
! polarized paraxial light beam carrying the orbital angular momentum 

suffers the 2D shift of the beam's centre of gravity relative the geo- 
. metric optic axis. The mutually orthogonal components of this shift 

qq | are expressed through the real and imaginary parts of the common 

O complex quantity. The features of the 2D vector, which describes the 

^ ■ shift, are analyzed. 
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I. Introduction 

In the communication [1], a new nonspecular effect, which accompanies the 
process of partial reflection of a paraxial light beam at a plane sharp interface 
of two isotropic transparent media, has been predicted. It was shown that, 
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if the incident beam carries the intrinsic orbital angular momentum (OAM), 
the centre of gravity of the reflected beam (CGRB) suffers a transverse shift 
(TS), i.e. the shift perpendicular to the plane of incidence. Bekshaev and 
Popov [2] considered the OAM-dependent TS using the method of transfor- 
mation of space-angle intensity moments during the interaction of the beam 
with a layered structure. Bliokh compared this effect with the phenomena 
accompanying transformation of the OAM in the course of propagation of a 
beam in the smoothly inhomogeneous media [3]. 

Unlike the previously known, spin-dependent TS (see, for instance, [4-9] 
and references therein), the OAM-dependent TS can take place, when the 
incident beam is s- or p-polarized. This effect appears because the electro- 
magnetic energy inside the beam is redistributed after reflection in such a 
way that the intensity increases at one side of the incidence plane and de- 
creases at the other [10,11]; in the regular region, the respective corrections 
are proportional to the ratio of the wavelength to the radius of the beam. 

The magnitude of the OAM-dependent TS is, if an angle of incidence 
is not close to the critical angle for total reflection or (for p-polarization) 
to the Brewster angle, rather small: it is of the order of the wavelength or 
less. But, at the present, the optical method has been developed, which 
permit to detect the positions of the objects with nanometer-scale precision 
(see, for instance, [12] and references therein). This method has been used 
by Dasgupta and Gupta [13] in order to define the OAM-dependent TSs of 
the CGRB in the wide region of the angles of incidence, the authors of [13] 
have reported a good agreement between their experimental results and the 
predicted values of the TSs. 

Recently, Okudo and Sasada have observed the redistribution of the elec- 
tromagnetic energy inside the beam carrying the OAM after its reflection 
from the transparent medium [11]. Their experiments were carried out in a 
situation, where the angle of incidence was close to the critical angle for total 
reflection. In this region, the magnitude of the effect increases significantly, 
so the energy redistribution becomes visible (see Figs. 3 and 4 in [11]). The 
authors of [11] have performed the respective numerical simulations, which 
turned out to be in a good agreement with the experimental results. The 
OAM-dependent TS was not calculated in [11], but this calculation was pos- 
sible [14]. 

Basing on the results of [11] and [13] one can conclude that the OAM- 
dependent TS of the CGRB is a quantity, which can be detected experimen- 
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tally, so, this effect can be used in order to investigate some features of the 
light and the matter. Thus, this effect is connected with the rotation motion 
of the electromagnetic energy inside the incident beam [15-17]; hence, if the 
OAM-dependent TS of the CGRB is detected, one can conclude that such a 
motion takes place. 

In this communication, we shall consider the OAM-dependent TS of the 
CGRB together with another shift, which is parallel to the plane of inci- 
dence and which corresponds to the well-known Goos-Hanchen shift (see, for 
instance, [5,8,9,18-21] and references therein); these shifts are of the same 
scale, moreover, they have a common origin. By this, we will assume that 
the reflecting medium may be lossy. 

II. Geometry of reflection 

We shall consider the reflection of a monochromatic light beam at a plane in- 
terface of two semi-infinite isotropic, nondispersive, and nonmagnetic media. 
The scheme of the process is shown in Fig. 1. The position of the interface is 
defined by equation N ■ p = 0, where p is the 3D radius vector, and N is the 
unit normal to the interface directed from the first (incident) medium to the 
second (reflecting) one. The dielectric constants of the media in the upper 
half-space (N • p < 0) and the lower half-space (N • p > 0) will be denoted 
by e x and e 2 , respectively. The relative dielectric constant e = e 2 /ei- The 
first medium is assumed to be transparent, while the second medium may 
be lossy, i.e. we assume that e\ is positive, while e 2 and, as a consequence, e 
may be complex, so, in general case, e = e' + ie". Here and further on, the 
prime and double prime stand for the real and imaginary parts of a quantity. 

Throughout this paper, we shall use the superscripts i and r for indicat- 
ing the quantities characteristic of the incident and of the reflected beams, 
respectively; the superscript a will be used in order to designate an arbitrary 
(incident or reflected) beam, so, a=i or r. We shall employ two coordinate 
systems connected with the incident and the reflected beams, these systems 
will be termed the i- and r-systems. 

In the i-system, the axis is assumed to coincide with the incident 
beam's axis. This axis and the unit vector N define the position of the beam's 
plane of incidence in 3D space; the angle 6 between the unit vectors N and 
zW is the beam's angle of incidence: 9 = arccos(N • z^). The coordinate 
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origin O is taken to be the point of intersection of the incident beam axis 
with the interface. In the r-system, the axis is assumed to coincide with 
the geometric optic (GO) axis of the reflected beam. The latter is defined 
as a ray, which is intersected by the interface at the coordinate origin, and 
whose direction is characterized by the unit vector z^ r \ which is related to 
zW through the Snell's law and is written as follows: = z^ — 2(N-z^)N. 
The beam's angle of reflection 9^ = arccos(N • z^) = ix — 9. 

The y axis in every system is defined to be perpendicular to the plane of 
incidence, it is characterized by the unit vector y = N x z^/|N x zW|; this 
axis (the transverse one) is common for both systems. The direction of the 
axis of the a-system, which lies in the plane of incidence, is defined by 
the unit vector x^ = y x z^ a \ 

In the a- system, the 3D radius-vector is represented as follows: 

p = z WzW + u (°> , (1) 

where and are the 2D planar radius vectors lying in the planes, which 
are perpendicular to the axis of the incident beam and to the GO axis of the 
reflected beam, respectively, 

u (a) = x (a) ± (a) + y ^ _ ( 2 ) 



III. Incident and reflected fields 

Let us assume that the incident beam is paraxial, that it carries the well- 
defined OAM ([22,23]), and that it is mainly s- or p-polarized. Due to first 
assumption, the relation 

\/nb < 1 , (3) 

must be fulfilled, where A is the wavelength of light in the first medium, and 
b is the beam radius. 

The electrical field vector of the incident beam under consideration can 
be written, in its own coordinate system, as follows [24]: 



E<?(P) ^ 



e« - i A*<0 f e« d 



2tt V a duW 



FW(p) + c.c, (4) 



where the subscript a denotes the polarization of the incident beam (a = s 
or p), while is the polarization vector: = y, and = xW. The 
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function F^\p) in the right-hand side of Eq. (4) can be represented as a 
product of two terms: 

F®(p) =/(«« expH/^W) , (5) 

where / is the azimuthal index (7 = 0, ±1,±2...), and y?W is the azimuth 
obtained from the relation: tan</?W = y/x®. The term f(u^\z^) does 
not depend on <£>W, for instance, it may be given by the Laguerre-Gaussian 
function or, like in [11], by a superposition of such functions. For the forth- 
coming analysis , the specification of f{vS l \z^) is not necessary. The time 
dependence of EW is suppressed. 

The electrical field vector of the reflected beam 'E^(p) is calculated in the 
standard way (see, for instance, [6,10,18]). First, the 2D Fourier transform 
of the function F (yt \p) with respect to the coordinates and y is to be 
performed. Then the vector E^(p) will be represented as a superposition of 
the plane waves. Next, the Snell and Fresnel laws should be applied to the 
particular plane waves. After that, the reverse 2D Fourier transform should 
be performed. 

Let us denote the angle of incidence of the particular plane wave inside 
the incident beam by 9^ p \ For the actual waves constituting a paraxial 
beam, \9^ -0| ~ A/fc<C 1. So, one can expand the particular Fresnel field 
reflection coefficients in the truncated power series with respect to 9^ — 9. 
When the second medium is transparent, and e' < 1, such an expansion can 
be substantiated, if 9 is not too close to the critical angle for total reflection 
given by 

9 ( - c) = arcsin V? , (6) 

i.e. if the condition 

\9-9 ic) \ > X/b, if e' < 1, and e" -> -0 , (7) 

is fulfilled. Again, in the case of p-polarization, the expansion is substanti- 
ated, if 9 is not too close to the Brewster angle given by 

9^ = arctanv 7 ^ , (8) 

i.e. if the condition 

\9 - 9 {B) | > X/b, if a = p, and e" -> -0 , (9) 
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is fulfilled. But, when 

|e"| > A/6 , (fO) 

the series expansion is justified at any 9. 

Let us restrict ourselves with the near-field region (in this case the possible 
angular shifts can be excluded from consideration), i.e. let us assume that 

Z^ < 7lb 2 /X . (11) 

If the restrictions (11), and (7), (9) or (10) are fulfilled, it is possible to retain 
only the zero-order and first-order terms in the power series expansion. 

Performing the operations, which are mentioned in three by last para- 
graph, one obtains the following expression for the electrical field vector of 
the reflected beam written in its coordinate system: 



E<;>(p) ~ rJ0) 



(12) 

where = y, = x^, and tan</?^ = y/x^'\ The quantity Q a (9) in the 
right-hand side of Eq. (12) is as follows: 



f(u {r) ,z (r) )exp(il<p {r) )+c.c. , 



Qa(e) = -~(\nr a {ff)) 



(13) 



where r a (i9) is the familiar Fresnel's field reflection coefficient of the a- 
polarized plane wave, which is incident at an angle i? [25]: 

. . sin(#-r) . . tan(#-r) 

r ' w = ~^IT7, ' r ' w = TM¥T7,- (14) 

and where r is the angle of refraction given by the relation y/e sin r = sin t?, 
in a general case r is complex. The dependence of the vector E^(/o) on 9 is 
not explicitly written down. 

Basing on Eqs. (13) and (14), the calculations of the factor Q s (9) and 
Q p {9) are straightforward, the results are as follows [10]: 
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and 

QpM = 1 Q 2 8 } 6) 2n ■ (16) 

e 1 sin 9 — cos 2 9 

Generally, Q a {9) is a complex quantity, its real and imaginary parts are 
caused by the dependence on 9 of the amplitude and the phase, respectively, 
of the Fresnel's field reflection coefficient r a (9). The expressions (15) and 
(16) are in accordance with the respective expression obtained in [21] (see 
Eq. (29)). 

The magnetic field vector of the reflected beam H^(p) is obtained by 
means of the Faraday's law. 



IV. 2D shift of the CGRB 

A. Definition and results of calculations . Let us select the plane, which 
is perpendicular to the GO axis of the reflected beam (the output plane in 
Fig. 1). The position of the CGRB on this plane (see Fig. 2) is given by the 
2D vector 

U«(0) = -4y / u««;M(pW r) , (17) 

where w£\p) is the electromagnetic energy density inside the reflected beam 

1 



«£>(p) 8 , 



ei(ESr>(p)) +(HW(p 



(18) 



and 

Wi r) = J «£>(p)duM . (19) 

As the GO axis intersects the output plane at the point u( r ) = 0, the 
vector XJ a (9), which is shown in Fig. 2, describes the 2D shift of the CGRB 
relative to this axis. 

Let impose the following restriction on the axial coordinate of the output 
plane: 

z {r) > b/ sin(2#) . (20) 

If the restriction (20) is fulfilled, the reflected field on the output plane is 
not overlapped with the incident field. Again, w^\u^) is negligible at the 
line of intersection of the output plane and the interface; in this case, the 
lower limit of integration over in Eqs. (17) and (19) can be substituted 
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by — oo. Once this substitution is made and the expression (12) for E^(p) 
and respective expression for H^(p) obtained by means of the Faraday's 
law are used, the calculation of U a (0)( r ) is straightforward. In the first-order 
approximation with respect to \/b, one obtains 

U a (0)=X a (0)*W+y a (0)y, (21) 

where 

*a(P) = QM~ , (22) 

TT 

and 

Y a (0) = IQ'M- ■ (23) 

TT 

B. Analysis. As follows from Eq. (21), the shift of the CGRB relative 
to the GO axis is, in a general case, described by the 2D vector \J a (9) on 
the output plane. Its a^ r ) -component (X a (9)) does not depend on I, while 
the y-component (Y a (6)) is proportional to I. Apart from the factor /, the 
functions X a {9) and Y a (8) are represented by the imaginary and real parts 
of the same complex quantity, it is in spite of the fact that the mechanisms 
of the appearances of the orthogonal components of the 2D shift are different 
[10,26]. Figs. 3(a,b) illustrate dependence of the functions Q' a and Q" a on 9 
for different values of e' and e" . 

When the second medium is transparent, the vector XJ a (9) is parallel 
either to the axis (in the total-reflection regime) or to the y axis (in the 
partial-reflection regime). In this case, the expression (21) reduces to the one 
of the previously obtained expressions, either for the Goos-Hanchen shift [27] 
or for the OAM-dependent TS of the CGRB [1]. When the second medium 
is lossy, the vector U a (9) is inclined to the coordinate axes, as is shown in 
Fig. 2. By this, the direction of U a (9) depends on the angle of incidence, 
and the vector U a (0) rotates around the coordinate origin, i.e. around the 
CO axis, when 9 changes. 

The direction of U Q (0) is defined by its azimuth T a (9), which is given by 
the equation: tanT a (9) = Y a {9) / X a {9). Substituting in its right-hand side 
Eqs. (22) and (23) one obtains the relation between the azimuth T a {9) and 
the argument of the complex quantity Q a (9), which is defined by j a (9) = 
(kiQ a (0))"); this relation looks as follows : 

tanT a (0) = Z/tan7 a (0) . (24) 
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In particular, if I — ±1, 



r«(fl) = /(|-7«W)- (25) 



As for the arguments 7 S (6 I ) and 7 P (0), they are obtained from the expressions 
(15) and (16): 

1 fe'-sm 2 9\ 
ls(v) = -arccotan — - 7r , (26) 



and 



/ e' — lei 2 cot 2 ^\ . „. 

7 P (0) = ls{0) ~ arccotan L±- . (27) 



Basing on Eqs. (24)- (27) one can make the following conclusions about 
the orientation of the vectors U s (9) and U p (9). As follows from Eq. (26), 
— 7r/2 > 7s(0) > — 7r at any 9, i.e. Qs(#) lies in the third quadrant of the 
complex plane; as a consequence, the vector U s (9) (see Eq. (24))lies in the 
third quadrant of the output plane (X s (9) < 0,Y S (9) < 0), if I > 0, and in 
the second quadrant (X s (9) < 0,Y S (9) > 0), if I < 0. The positions of Q p {9) 
on the complex plane and of the vector \J P (9) on the output plane, unlike 
of Q s {9) and U s (9), are not contained within one quadrant, when 9 changes 
from to 7r/2. At small 9 the direction of the vector U p (9) is approximately 
opposite with respect to the direction of the vector U s (9). As 9 increases, 
the sign of Y p {9) changes once independent of the value of e. The sign of 
X p {9) changes, if the condition 

2e' > |e| 2 (28) 

is fulfilled. The numerical analysis of dependence of X on 9 has been recently 
performed by 

The angles, at which these changes take place, are obtained from the 
following equation: 

|e| 2 cot 4 # + 2(l + e'|er 2 )shi 2 #-3 = . (29) 

If e" — > —0 and e' > 0, the second term in the right-hand side of Eq. (27) 
have the jump discontinuity at the Brewster angle 9^ B \ where 7 P changes 
abruptly by n. If in addition e' < 1, the (7r/2)-jump of 7 S and 7 P takes place 
at the critical angle for total reflection 9^ c \ Figs. 4(a,b) illustrate dependence 
of the functions T p and T s on 9 for different values of e and /. 
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The scale of the 2D shift is X/tt, but, if \e"\ <C 1, its magnitude increases 
significantly in the vicinity of the angle 9^ and, for p-polarization, in the 
vicinity of the angle 0( B \ these angles are given by Eqs. (6) and (8), re- 
spectively. Fig. 5 illustrates the behavior of the vector U s (9) in the former 
domain, the behavior of the vector U p (0) in this domain is similar (see fur- 
ther Eq. (30)). It is seen that, when 1^0, the effect of rotation of the vector 
U a (9) around the GO axis as well as change of its length in the process of 
the change of 9 is pronounced. Again, if \e"\ <C 1, and \9 — 9^\ <C 1, the po- 
sition of the vector XJ a (9) is sensitive to small changes of the values of e' and 
e" (see Figs. 6 and 7). If the above-mentioned conditions are fulfilled, the 
approximate expressions of the functions Q' a {9) and Qa(9) look as follows: 

«"W = WW ="( 5^ ] l 1 -^^) 1 ' 2 ■ (30) 




where A — (9 — # {c) ) sin(26 l(c) )/ |e"| , in (30) and further on the signs minus 
and plus stand for Q' a (9) and Q'^(9), respectively. Notice that, relative to the 
point 9 = 9^ c \ i.e. to the point A = 0, the functions Q' a (9) and given 
by Eq. (30) are mutually symmetric. The magnitudes of these functions 
achieve the maximums at the symmetrical points A^. and A^,, respec- 
tively, given by = -A^ = l/y/3. By this, Q' s (A$ ax ) = Q' S '(AW J = 

— (3/4) 3 / 4 (eVje ,, j) 1/2 - 

The behavior of the vector XJ a (9) in the vicinity of 9^ is shown in Fig. 8. 
Like in the previous case, the effect of rotation of this vector around the GO 
axis with the change of 9 is pronounced in this domain; again, the position 
of U p (0) is sensitive to small changes of the values of e' and e" (these effects 
are not demonstrated). When |e"| < 1, and \9 - 0< B > | < 1, 

Q w*M s -^ +i w?hT ' (31) 

while |<3s(#)| is small in comparison with \Q p (9)\, namely, Q' s (9) ~ — l/y/e*, 
and \Q't(9)\ « \Q' S (9)\- 

The numerical analysis of dependence of Goos-Hanchen shift, i.e. of the 
quantity X, on 9 has been recently performed by Lai et al. The authors of 
[21] have paid special attention to the behavior of this shift in the vicinity 
of the Brewster dip. Eq. (31) shows that, if \e"\ 1, the function X p (9) is 
approximated in this region by the Lorentz curve. 
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C. Relative p-s shift. Detection of the shift of the CGRB of the s- 
or p-polarized beam requires a careful determination of the position of the 
GO axis. The latter procedure can meet the difficulties. One can escape 
these difficulties by operating with the difference of the shifts of two beams 
having the orthogonal polarizations, such a scheme has been realized in the 
experiments by Dasgupta and Gupta [13]. In [13], the initial beam was 
randomly polarized, and the s- or p-polarization were selected by means of 
rotation of a linear polarizer. The other elements of the set-up were kept 
motionless. 

/,From Eqs. (17) and (19), the difference between the centres of gravity 
of the respective reflected beams 

\J ps (9) ee \J P (9) - XJ S (9) = ! f u W(?W(uW)duW, (32) 

where GW(uW) is the cross-correlation function, 

GW(u«) = J ic«(u (r) + v( r Vi r) (v (r) )rfv (r) , (33) 

and where is, like u^, the 2D radius vector lying in the plane perpen- 
dicular to e^. The vector U ps (9), unlike the vectors U s (6>) and U p (6>), does 
not depend on the position of the GO axis. 

Let us consider the case, when both incident beams have the common 
axes and, as a consequence, both reflected beams have the common GO 
axes. Such a situation takes place, for instance, if, within the scheme of the 
experiment [13], the rotation of the polarizer do not affect the position of the 
axis of the incident beam. In this case 

US(fl) = (W (r) + W)^ (34) 

where 

Q P s(0)=Q P (0)-Qs(0) . (35) 

Generalization of the expression (34) for the case, when the axes of the 
mutually orthogonally polarized incident beams are parallel, is straightfor- 
ward. 
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V. Conclusions 



When the s- or p-polarized paraxial light beam, which carries the OAM, is 
reflected from the lossy medium, two mutually orthogonal shifts of the CGRB 
to scale A appear; so, in the near-field region, the total shift is described by 
the 2D vector XJ s (0) or U p (6>) on the output plane. The vector U a (9), with 
a = s or p, is given by Eq. (21) and its orthogonal components are given by 
Eqs. (22) and (23). The second medium is transparent, these expressions are 
correct, if conditions (7) and (9) are fulfilled. The second medium is lossy, 
and condition (10) for \e"\ is fulfilled, the above-mentioned expressions are 
correct at any 9. 

The vectors \J S (9) and U p (9) are defined relative the GO axis of the 
reflected beam; the difference between them (the vector U ps (9) = U p (9) — 
XJ S (9)) does not depend on the position of the GO axis, it is expressed through 
the cross-correlation function of two energy distributions across the output 
plane (Eqs. (32) and (33)). 

The component of the vector XJ a (9), which is parallel to the plane of 
incidence (X a (9)) corresponds to the Goos-Hanchen shift, it does not depend 
on the azimuthal index /. The transverse component Y a (9), which represents 
the OAM-dependent TS of the CGRB, is proportional to /. Apart of the 
factor /, the components X a {9) and Y a (9) are given by the imaginary and 
real parts of the common quantity. 

Redistribution of the electromagnetic energy inside the reflected beam, 
which entails the OAM-dependent TS of the CGRB, is connected with the 
rotation energy motion inside the incident beam [15-17]; hence, if a deflection 
of the vector U Q (#) from the x^> axis is detected, one can conclude that such 
a motion takes place. 

In the case of s-polarization the signs of both components of 2D shift do 
not depend on 9: the vector XJ S (9) lies in the third quadrant of the output 
plane, if / > 0, and in the second quadrant, if / < 0. On the contrary, the 
positions of the vector ^J p (9) on the output plane is not contained within 
one quadrant when 9 changes from up to n/2. The azimuths of the vectors 
U s (0) and U p (9) are given by Eqs. (24)- (27). At small 9 the directions 
of these vectors are approximately opposite . As 9 increases, the sign of 
Y p (9) changes one time independent of the value of e, while the sign of X p (9) 
changes, if the condition (28) is fulfilled. Generally, when 9 changes, the 
rotation of the vector XJ a (9) around the GO axis as well as the change of 
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its length takes place; these effects are most pronounced if \e"\ <C 1 and 
9 is close to 9^ or, for a — p, to 9^ B \ In these domains, the quantity 
Q a (9) is approximated by expression (30) or (31). In the vicinity of 9 (c \ the 
dependence of the vector XJ S (9) on 9 is shown in Fig. 5, the behaviors of 
the vectors U p {9) and XJ ps (9) are, due to the left relation in (30), similar to 
behavior of XJ S (9). In the vicinity of 9^ B \ the dependence of the vector U p (#) 
on 9 is shown in Fig. 8. 

In the above-mentioned domains, the position of the vector XJ a (9) is sen- 
sitive to small changes of the values of e' and e" (see Figs. 6 and 7). This 
may be used as a tool for the investigation of the small changes of the di- 
electric constants of the media. Again, in view of this fact, one can await 
that, if the reflecting medium is taken to be nonlinear, the dependence of the 
position of the vector XJ a (9) on the light intensity should be pronounced in 
the above-mentioned domains. 
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FIGURE CAPTIONS. 

Fig. 1. Geometry of reflection. 

Fig. 2. Position of the vector U on the output plane, the azimuth T is 
shown, the subscripts a are suppressed in this figure. The numerals 1, 2, 3, 
and 4 denote the quadrant' numbers. 

Fig. 3(a). Q' a (thin curves) and Q" a (thick curves) versus 9. e' = 2; solid 
lines correspond to e" = —0.001, and dashed lines to e" = —0.5. 

Fig. 3(b). Same as in Fig. 3(a) but with e' = 0.5 

Fig. 4(a). T a versus 9. e' = 2. Solid curves for I = 1, e" = -0.001, 
dashed curves for I = 1, e" = —0.5, and dotted curves for I = 5, e" = —0.5. 
Fig. 4(b). Same as in Fig. 4(a) but with e' = 0.5 

Fig. 5. Dependence of the components X s and Y s on 9 in the vicinity 
of 9 (c K e' = 0.5 (0< c > = 45°), and e" = -0.001. The thin, middle, and thick 
curves correspond to I — 0, 1, and 2, , respectively. 

Fig. 6. Dependence of the components X s and Y s on e' in the vicinity 
of e' = sin 2 #. I — 1, 6 — 45°. The thick and thin curves correspond to 
e" = —0.001 and e" = —0.01, respectively. 

Fig. 7. Dependence of the components X s and Y s on e" in the situation, 
when 9 is close to 9^. I — 1, e' — 0.5. The thin and thick curves correspond 
to 9 = 44.9° = # (c) - 0.1° and 9 = 45.1° = # {c) + 0.1°, respectively. 

Fig. 8. Dependence of the components X p and Y p on 9 in the vicinity of 
9 (B \ e' = 2.0 (9^ = 54.7°), and e" = -0.001. The thin and thick curves 
correspond to / = 1, and 2, respectively. 
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